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O ■ Abstract 



We describe a pure state of four qubits whose single-qubit density 
matrices are all maximally mixed and whose average entanglement as 
a system of two pairs of qubits appears to be maximal. 

a 

1 Introduction 

> ■ 

• i-H , 

^ ■ In a system of two qubits the state 

i 

|C 2 ) = - 7 =(|00) + |11)) (1.1) 

is, on all counts, the most entangled of all pure states. It gives the great- 
est violation of Bell inequalities, it has the largest entropy of entanglement, 
and its one-party reduced states are both maximally mixed. All of these 
properties determine it uniquely up to local unitary transformations. 
A pure state of three qubits with similar properties is the GHZ state 

|C 8 > = -^(|000) + |111». (1.2) 

This state has the maximum value of pure 3-party entanglement, as measured 
by Wootters's 3-tangle [|, ||, and its one-particle reduced density matrices 

*See ref. 0. 



1 



are all maximally mixed. Like the two-qubit state IC2), it is characterised 
uniquely, up to local unitary transformations, by the latter property ||. 
The obvious n-party generalisation is the "Schrodinger cat" state 

|C„> = -^(|00...0) + |11.-.1». (1.3) 

Like \Cz) and IC3), this state has the property that its one-party reduced 
states are all maximally mixed. On the strength of this, \C n ) is sometimes 
called the "maximally entangled" pure state of n qubits. For n > 3, however, 
not all states with this property are locally equivalent, and it is not clear that 
\C n ) is really the most entangled of them. Here we examine the case n — 4, 
show that there are 4-qubit states which are more entangled than IC4), and 
attempt to find the most entangled among them. 

A four-qubit system can be regarded, in three different ways, as a system 
of two pairs of qubits, and one can ask how entangled are these pairs. In the 
state IC4) this entanglement is not maximal: each pair X, Y of the four qubits 
A,B,C,D is not in the maximally mixed state but exhibits correlations, all 
two-qubit density matrices being 

Pxy = ±(|00>(00| + |11>(11|). (1.4) 

But there do exist pure states of four qubits in which four of the six two-qubit 
reduced density matrices are maximally mixed. (We note that these density 
matrices come in pairs: if the reduced state of one pair is maximally mixed, 
so is that of the complementary pair.) For example, the state 

= I (|0000) + |0111) + |1001) + |1110» (1.5) 

has two-qubit density matrices 

Pab = Pac — Pbd = Pcd = jl, 

p^ = !(| + +X+ + | + |--X--|), 

Psc = I (|00)(00| + |11)(11|) , 



where 

|±) = -^(|0>±|1». (1.6) 

The one-qubit reduced density matrices of are all maximally mixed. But 
|\&) is not locally equivalent to the cat state IC4), for it has different values 
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from that state of the entanglement entropies Eab, Eac, Ebd, Ecd, which 
are invariants under local unitary transformations. (We write 

E X y = -tr(p XY log 2 p X y), Pxy = tr zw |*)(*| (1.7) 

where {W, X,Y, Z} is a permutation of {A,B,C,D}). In fact QJ, the state 
|\&) cannot be asymptotically reversibly converted into any collection of the 
states \C n ) for n = 2,3,4. 

It is natural to wonder whether there is a still more entangled state in 
which the reduced density matrices of all six pairs of qubits are maximally 
entangled.Q In Section 2 we will show that there is no such state of four 
qubits. The situation changes, however, if particles with larger state spaces 
are considered, and we will show that a system of four four-state particles 
does have pure states with this property. 

In Section 3 we ask what is the maximum two-pair entanglement possible 
for a system of four qubits. Taking as a measure the average entanglement 
entropy (E 2 ) = ^(Eab + Eac + Ead), we exhibit a pure state |M 4 ) with a 
greater value of this quantity than the state (|1.5f) ; we show that (E2) has a 
stationary value at |M 4 ), and present evidence that this is a global maximum. 

2 Non-existence of maximal entanglement 
between all pairs of pairs 

Theorem 1. There is no pure state of four qubits whose two-qubit density 
matrices are all multiples of the identity. 

Proof. Write the four-qubit state |\&) as 

i*>= E tl3kl \m\m) 

i,j,k,l=0,l 

The tensor t l ^ kl can be regarded as a 4 x 4 matrix in three different ways: 

and the requirement of maximal entanglement is that U\, U2 and U3 should 
all be unitary matrices. 

lr This question has been studied by Gisin and Bechmann-Pasquinucci Q for the general 
case of n qubits, but under the restriction that the states are symmetric under permu- 
tations of the qubits. This is an unnatural requirement in this context, since it is not 
invariant under local unitary transformations. 
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We show first that by local unitary transformations we can arrange that 
the coordinates of satisfy 



,1000 



,0100 



,0010 



,0001 



0. 



To do this, we find normalised states \a), \/3), I7), \S) which maximise N = 
I j ck) j^) j <5) j 2 . Such states certainly exist, since N(a, (3, 7, 5) is a contin- 
uous function on the compact space S 3 x S 3 x S 3 x S 3 . Change basis in 
each of the single-qubit spaces so that |a)|/3)|7)|<5) becomes the basis state 
|0000); then N = |t 0000 | 2 . But if t 1000 were non-zero we could change basis 
for states of the first qubit so as to increase |t 0000 |; since we have maximised 
it, we must have t 1000 = 0. Similarly t 0100 = t 0010 = t 0001 = 0. 
The matrices Ui, U 2 , U 3 which have to be unitary are now 



A 1 



V 



0000 











t owi 


£0110 





£1001 


£1010 


1100 


t noi 


£1110 



£0011\ 
£0111 
£1011 
£1111 



V 



0000 











£0011 


£0110 





£1001 


£1100 


1010 


£1011 


£1110 



£0101\ 
£0111 
£1101 
£1111 



V 



0000 








£0110 





£0011 


£0101 


£0111 





£1010 


£1100 


£1110 


1001 


£1011 


£1101 


£1111 



Hence the following three conditions must hold: 

1. t°on = or t 0111 = t 1011 = 0; 

2. t 0101 = or t 0111 = t 1101 = 0; 

3. t 0110 = or t 0111 =t 1110 = 0. 

Following the branching consequences of these, and requiring that the ma- 
trices are unitary, leads to a conclusion of the form that three 2x2 matrices 



a b 
c d 



b e 
f c 



a e 
f d 



must all be unitary. In this situation all the matrix elements must be non- 
zero: if, for example, a = 0, then d = and b, c, e, f would all have unit 
modulus, which is impossible if the second matrix is to be unitary. Now 
unitarity gives 



bd 



a = — - 



ed 

T 



so 



But 



bf = ec. 



bf + ec = 0, 
so bf = ec = 0, which is impossible. 



□ 
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Four-party states with this property do exist, however, if the individual 
parties have more independent states. For example, suppose each party has a 
4-dimensional state space; then the (non-normalised) state with coordinates 

f 1 if i = j = k = I or k, I) is an even permutation of (1, 2, 3, 4) 
1 otherwise 

has every two-qubit density matrix equal to a multiple of the identity. 

We note that the method introduced in this proof can be used to define a 
canonical form for any multipartite pure state, in which a maximal number of 
coordinates are zero. This can be regarded as a generalisation of the Schmidt 
decomposition of a bipartite state. Details can be found in 0. 



3 A Maximally Entangled Four-qubit State? 

Given that a four-qubit state cannot have maximal entropy of entanglement 
for every two-qubit subset, we now ask what is the greatest possible average 
for such entropies, i.e. we seek to maximise 

(E 2 ) = \{E AB + E AC + E AD + E BC + E BD + E C d) 
= \{Eab + Eac + Ead)- 

The second equality holds because complementary pairs have equal entropy. 

We have not been able to solve this problem analytically. We will adopt a 
heuristic approach, using an (indefensible) analogy with the two-qubit system 
to obtain a candidate maximally entangled state and then showing that (E 2 ) 
is indeed stationary at this state and appears to be maximal. 



Maximally entangled states of two qubits like ( PD are often referred to 
as "singlet" states, since an example of such a state is the state of two spin-| 
particles with zero total angular momentum. This description is not invari- 
ant under local unitary transformations, which need not preserve the total 
angular momentum. Nevertheless, let us take this as a hint in investigating 
four-qubit states. A system of four spin-1 particles has two independent sin- 
glet states. The most symmetric combination of these, in the sense that all 
of its two-qubit density matrices are unitarily equivalent, is locally equivalent 
to 

|M 4 ) = —= [|0011) + |1100) + cjfllOlO) + |0101)) + a/(|1001) + |0110))1 , 
V6 

where to = e 2m / 3 , or to |M 4 ), the complex conjugate of |M 4 ) in the compu- 
tation basis. Note that this state is not symmetric between the qubits but 



5 



belongs to a two-dimensional representation of the permutation group, the 
other state in the representation being |M 4 ). 
The two-qubit density matrices of |M 4 ) are 

Pab = PAC = PAD = \ [|00)(00| + |11)(11| + |$+)($+|] + (3.1) 

where 

|*±> = i(|10>±|01» 

(in angular momentum terms, the two-qubit reduced states are equal mix- 
tures of a singlet and a maximally mixed triplet). Hence the entanglement 
entropies are 

Eab = E AC = E AD = 1 + | log 2 3. 
Comparing with the cat state IC4), for which 

Eab = Eac = Ead = 1, 
and the state of (1.5Q, for which 

Eab = Eac = 2, Ead = 1, 

we see that \M^) has a greater value of (E 2 ) than either. 

We will now show that the function (E 2 ) is stationary at |M 4 ). To simplify 
the calculation, we consider the functions E X y defined by ( |1.7| ) for all state 
vectors \^f), though these functions coincide with the entropies only when 
|\&) is normalised. Suppose changes by \S^f). Because the trace makes 
the operators behave as if they commuted, the consequent change in Exy is 
given by 

5 E XY = - tr [p X v log 2 (ep X y)] • 
At = VE\M 4 ) we find 

5E AB = Re [-(*!<**) log 2 (3e 2 ) + log 2 3] , 

5E AC = Re [-(*|«y*> log 2 (3e 2 ) + u<p\6V) log 2 3] , 
5 E AD = Re [-(^\8^f) log 2 (3e 2 ) + o; 2 (¥|5*) log 2 3] , 

where \^f) is the complex conjugate of |\&) in the computation basis. Hence 
the gradient of the average (E 2 ) in the real 32-dimensional space of all state 
vectors is in the direction of |^): 

5(E 2 ) = -log 2 (3e 2 )Re(^|<^). 
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Thus (E2) is stationary at \^) for variations on the sphere of vectors with the 
same norm as But each Exy, as given by (|1.7|) , is linearly related to the 
true entropy of entanglement, so the average two-party entropy is stationary 
at |M 4 ). 

We have searched numerically for states which maximise (E 2 ) starting 
from several arbitrarily chosen states. All the states we have obtained in this 
manner are locally equivalent to IM4) or IM4). 



4 Robustness of the entanglement 

One of the criteria for maximal entanglement suggested in J3J is that the 
entanglement should be maximally fragile: a measurement on any one of the 
qubits destroys the entanglement between the remaining qubits. This is true 
of the cat state (|1.3|) if the measurement projects onto the computation basis 



{|0),|1)}. Projection onto the basis {|+), |— )} defined in ( |L6D , however, 
leaves the remaining qubits in a cat state. 

The state constructed in the previous section behaves in the opposite way 
to the cat state: measurement of one qubit leaves the other three qubits in an 
entangled state, and the amount of entanglement afterwards is independent 
of what measurement is performed. To be precise, projection of one qubit 
onto a computation basis state leaves the remaining qubits in one of the 
entangled states 

^=(|011)+u;|101)+^|110)) (4.1) 

or —= (|100) +w|010) +w 2 |001)) . (4.2) 
v3 

These states, which are clearly locally equivalent, have recently been iden- 
tified as having maximal average two-qubit entanglement in a certain sense 
||. Projection onto any other state of the first qubit leaves the other three 
qubits in another state which is locally equivalent to the above. This follows 
from the fact that the state |M 4 ) is an SU(2) singlet, so that the unitary 
transformation from the computation basis to any other basis of one qubit 
can be undone by performing the same transformation on the other three 
qubits. 

Thus the entanglement of the state I-M4) is robust. Any carelessness by 
the holder of any one of the qubits, resulting in an uncontrolled decoherence 
of that qubit, does not completely destroy the entanglement of the remaining 
qubits, and always leaves them with the same amount of entanglement. 
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